DESY 95-077 
April 1995 



ISSN 0418-9833 



A Parton Model for Diffractive Processes 
in Deep Inelastic Scattering 

W. Buchmiiller and A. Hebecker 

Deutsches Elektronen- Synchrotron DESY, 22603 Hamburg, Germany 



Abstract 

We demonstrate that the global properties of the "rapidity gap" events, observed 
at HERA, can be understood based on electron-gluon scattering as the underlying 
partonic process. Using the measured inclusive structure function F2 to determine 
the parameters of the parton model, the diffractive structure function F2 is pre- 
dicted. The ratio of diffractive and inclusive cross sections, Rd = OBjoind ~ 1/9, 
is determined by the probability of the produced quark-antiquark pair to evolve 
into a colour singlet state. This colour singlet cluster may fragment into hadrons 
independently of the proton remnant, yielding the observed gap in rapidity. 



In the "rapidity gap" events, observed in deep inelastic scattering at HERA 
the detected hadronic final state has small invariant mass, and it is separated by a gap in 
rapidity from the proton beam direction. The absence of a hadronic energy fiow between 
proton remnant and current fragment suggests that in the scattering process a colour 
neutral part of the proton is stripped off which fragments independently of the proton. 
In analogy to hadronic processes of similar kind the "rapidity gap" events are also called 
"diffractive" events. 

The measured "diffractive" cross section is not suppressed at large values of the 
momentum transfer relative to the inclusive cross section. Such a "leading twist" 
behaviour is usually regarded as evidence for scattering on point-like constituents. This 
interpretation, however, appears to be in conflict with the fact that quarks and gluons, 
the constituents of the proton, carry colour. Hence, one expects the formation of jets in 
the flnal state with a hadronic energy flow familiar from ordinary deep inelastic scattering 
events without a rapidity gap. 

In the following we shall demonstrate that this puzzle can be resolved by taking non- 
perturbative fragmentation effects into account. We start from the production of a quark- 
antiquark pair in electron-gluon scattering as basic partonic process. Immediately after 
their production quark and antiquark propagate in the colour fleld of the proton. With 
a certain probability, approximated by a statistical weight factor, the quark-antiquark 
pair evolves into a colour singlet parton cluster which can escape from the proton and 
fragment independently into hadrons. The cross section for diffractive events can then be 
calculated in terms of the statistical weight factor, the quark-antiquark production cross 
section and the gluon density which we determine in terms of the inclusive structure 
function F2. Our approach is related to previous work on "aligned jet models" @]-[^ and 
"wee parton lumps" in deep inelastic scattering If. 

Quark-antiquark pairs are produced in electron-gluon scattering. The relevant kine- 
matic variables are (cf. flg. 1) 



which characterize inclusive deep inelastic scattering, and the invariant mass M of the 
quark-antiquark pair. 
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M' = {q + pgf. 
one has. 



(2) 



(3) 



2 



The differential cross section for the inclusive production of quark-antiquark pairs is 
given by [|] 

da{ep e{qq)X) a /..((^ , f , \ \ fA\ 

Here g{^) is the gluon density, and the cross section aT{L) is obtained by integrating the 
differential parton cross section over the momentum transfer t = {q — l'Y (cf- fig- 1) from 

tmin ~ Q / /3 to trraax ~ f^'' 

"^<« = I — ■ 

The differential parton cross sections (see fig. 1 and the crossed process) read 
dar iraas 9^9/,^/. ^^ u t , ^9, 

9 



dt 



where the sum extends over all quarks whose mass is small compared to Q. 

From eqs. @) and (|]) one obtains the contribution to the inclusive structure function 

A(^)F2(x,Q^)=xfiEe;rf^?(0f(/?^ + (l-/?r)ln4^-2 + 6/?(l-/3)) . (8) 

This is the result in the "massive gluon" scheme The virtuality —vn^g of the gluon 
regularizes the coUinear divergence at t = 0. In the more familiar minimal subtraction 
scheme one obtains a similar expression, where nig is replaced by the subtraction scale 
/i, and the finite part is different. At one-loop order, the complete expression for the 
structure function F2 reads, 

g^) = X 5: e\ C f (5,(0 + 5,(0) m - /5) + 0{a,)) + A(^)F2(x, Q^) . (9) 

Here Sg and Sg are quark density and antiquark density, respectively. 
At small values of ^ one has, 

5,(0 , 5,-(0 « 9(0 ■ (10) 
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For simplicity, we shall neglect in the following the quark contribution to F2. In this 
approximation we shall obtain a parameter free prediction for the diffractive structure 
function. Note, that Sq{^) is different from Sq{C,,Q'^), the quark density at scale Q^. 
Neglecting Sq{^) in eq. (D essentially amounts to calculating Sq{^,Q^) in terms of an 
"intrinsic" gluon density g{C,), corresponding to a scale rrig = 0{1 GeV^). A more com- 
plete analysis keeping the contribution from Sq{^) would require information about the 
longitudinal structure function Fi in order to predict the diffractive structure function. 

For the gluon density at small values of C, we use the usual parameterization, 

g{0 = A, r'-' , (11) 

where Ag is a constant. Inserting eq. (|ll|) into eq. (§) we can now evaluate the inclusive 
structure function F2, 

F2{x,Q') = ^J24^9ix)J^'dpp'(^iP' + il-Pr)\n^-2 + 6P{l-P)] . (12) 



A simple approximation for F2, valid at small x, is obtained by choosing x = as lower 
limit of integration, and by setting A = in the integrand. This yields (x <^ 1), 

F2{x, Q^) - E 4^9i^) il + In ^1 ■ (13) 

This expression provides a good description of the HI measurement of the structure 
function [|TT| for A = 0.23, nig = 1.0 GeV and AgasJ2q^l ~ 0.61. For these parameters 



(|T3p is very close to the phenomenological fit of the structure function F2 given in 
We have also evaluated the expression (^) without any approximation. Instead of the 



simple result (|T3|) one then obtains some cumbersome function of x and A. Comparison 
with the HI data yields again A = 0.23, whereas the value of rrig is now 0.8 GeV. The 



following discussion will be based on the simple analytical expression (13) 



We are now ready to calculate the diffractive structure function. The main idea is that 
the quark-antiquark pair, originally produced in a colour octet state, changes its colour 
randomly through further soft interactions with the colour field of the proton remnant. 
Hence, the quark-antiquark pair evolves into a parton cluster which separates from the 
proton remnant with some probability Pg in a colour octet state, and with probability 
Pi = 1 — Pg in a colour singlet state. In the first case, there is a colour flow between 
proton remnant and current fragment leading to the typical hadronic final state. In the 
latter case, however, the colour singlet final state may fragment independently of the 
proton remnant, yielding a gap in rapidity. For a sufficiently fast rotation of the colour 
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spin of quark and antiquark, the probabilities should simply be given by the statistical 
weight factor accounting for the possible states of the quark-antiquark pair, i.e., 

P,^l,Ps^l. (14) 
Similar ideas concerning the rotation of quarks in colour space have been discussed by 



Nachtmann and Reiter in connection with QCD-vacuum effects on hadron-hadron 
scattering. 

In analogy with the usual inclusive structure functions, the diffractive structure func- 
tions are defined as, 

'™'^^-y + ^Vf(-.Q^O-^^■i'(-.Q^o) . (15) 



The structure function is easily obtained from eq. (^. With x = and including 
the statistical weight factor (H), one obtains 

Fi'{x,Q\0 - l^T.<9{OF,''{f3,Q') , (16) 



where 

F2^(/?, Q2) = p f^p^ + (1 _ ;5)2) In ^ _ 2 + 6/3(1 -/?)). (17) 



"gf- 

Since the gluon density g{C,) and the mass scale nig have been determined by the fit to 
the inclusive structure function F2, the diffractive structure function is unambiguously 
predicted, including its normalization. 

Let us discuss the main properties of this result. An immediate consequence is the 



prediction of the ratio of diffractive and inclusive cross sections. From eqs. (^-(10) and 
(|16D one obtains, 

F,(x,Q^) -9- 

Note, that within the model described above, this ratio directly measures the probability 
of forming a colour singlet parton cluster in the scattering process. 



The form of the diffractive structure function (16) is identical with expressions ob- 
tained based on the idea of a "pomeron structure function" |12]-[l7]. The interpretation 
of the ingredients, however, is rather different. The "pomeron flux factor" is replaced 
by the density of gluons inside the proton, which factorizes. The "pomeron structure 
function" for partons with momentum fraction (3 inside the "pomeron" is identified as 
the differential distribution for the production of a quark-antiquark pair with invariant 
mass M2 = Q2(i _^)/^. 
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The function F^dSjQ"^) is plotted in fig. 2 for three different values of Q^. For 
intermediate values of f3 between 0.2 and 0.6 this function is rather flat. Approximating 
^FiP^Q"^) this interval by ^2^(0. 4, Q^), a comparison of eqs. and (p^) yields the 
scaling relation |p 

Fi^{x,Q\0^jm,Q') , (19) 

where D ~ 0.04, independent of Q^. Note, that the /5-spectrum shown in fig. 2 is rather 
sensitive to the infrared cutoff nig at small and large values of /3. 

The scaling relation ([T9|) provides a rather accurate description of recent measure- 
ments of the diffractive structure function by the HI collaboration. The experimental 
data are consistent with 

Ff(x,Q^O oc ln(Q2)^-, (20) 

where n = 1.19±0.06±0.07 0. This is in good agreement with our results eqs. {^) , (p^) 
with A = 0.23. 

Let us finally verify that the model described above predicts indeed the appearance 
of rapidity gaps. For simplicity, we shall consider the rapidity of the antiquark with 
momentum / in the 7*p-rest frame (cf. fig. 1). It is related to other kinematical variables 
by 

u + rn?l3 



77 = - In 
' 2 



e(l-/5) 



(21) 



From eq. (^) one obtains, 

dar) dt dan 



dydQ'^d^drj drj dydQ'^d^ dt 

dtl a y'^\(d(TT , daA y"^ daA , . 

= Tv-9^^'^^^[['-y-Y)[-dr + ^)-j^) ■ ^^'^ 

The total diffractive cross section for a maximum rapidity rimax can now be obtained 
by integrating over a specified kinematic domain, where the rapidity of the antiquark is 
larger than the rapidity of the quark. The reverse configurations yield the same contri- 
bution, resulting in a factor of 2. Hence, one obtains 

2 ,Q'' r t" ^i^^^ . (23) 



dr]max JqI Jyi Jx dydQ'^d^dr] 

where ^rnax = minjl, X + exp 277max}- Using the kinematic boundaries 0.03 < y < 0.7 and 
7.5 GeV^ < < 70 GeV^ one obtains the distribution shown in fig. 3. Note, that 
the approximation <^ has been used throughout the calculation. Above rjmax ~ 2 
the diffractive cross section is clearly negligible. Assuming that the final state hadrons 



6 



are produced in the rapidity interval spanned by the quark- antiquark pair, one obtains 
the observed rapidity gap. 

In summary, we have demonstrated that electron-gluon scattering can account for the 
global properties of the "rapidity gap" events observed at HERA provided the following 
two hypotheses concerning the formation of the final state are correct: First, the initially 
produced quark- antiquark pair evolves with a probability, given by a statistical weight 
factor, into a colour singlet parton cluster which fragments independently of the proton 
remnant; second, the rapidity range of the diffractive hadronic final state is essentially 
given by the rapidity interval spanned by the produced quark-antiquark pair. 

This simple picture appears to provide a rather accurate description of the observed 
diffractive events, including the total rate, the ,^-dcpcndcncc and the Q^-dcpcndcncc. This 
agreement with experimental data may appear fortuitous, since a number of theoretical 
issues still remain to be settled. These include the non-perturbative formation of the 
colour singlet cluster and the role of the infrared cutoff. Our results indicate, that an 
appropriate starting point for the evaluation of the inclusive as well as the diffractive 
structure functions at small values of x may be a semiclassical approach where "wee 
partons" are treated as a classical colour field. We expect that such a framework will 
lead to results very similar to the ones described above. 

We would like to thank G. Ingelman, H.-G. Kohrs, M. Kramer, M. Liischer and 
P. Zerwas for valuable discussions. 
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Figure captions 

Fig.l Quark-antiquark pair production in electron-gluon scattering. 

Fig. 2 Dependence of the diffractive structure function on /3 and Q^. 

Fig. 3 Distribution of the maximal rapidity in the 7*p-rest frame defined by the most 
forward quark or antiquark. 
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